(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     20914,        576]
NotebookOptionsPosition[     19858,        537]
NotebookOutlinePosition[     20246,        554]
CellTagsIndexPosition[     20203,        551]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"Clear", "[", 
  RowBox[{"GN", ",", "J", ",", "eta", ",", "psi"}], "]"}]], "Input",
 CellChangeTimes->{{3.719372911104762*^9, 3.7193729212845073`*^9}, {
  3.7193739044971285`*^9, 3.7193739073079395`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"GN", ":=", 
   RowBox[{".25", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"eta", "-", "1"}], ",", 
        RowBox[{"1", "-", "eta"}], ",", 
        RowBox[{"1", "+", "eta"}], ",", 
        RowBox[{
         RowBox[{"-", "eta"}], "-", "1"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"psi", "-", "1"}], ",", 
        RowBox[{
         RowBox[{"-", "psi"}], "-", "1"}], ",", 
        RowBox[{"1", "+", "psi"}], ",", 
        RowBox[{"1", "-", "psi"}]}], "}"}]}], "}"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.7193717611291685`*^9, 3.7193718112667513`*^9}, {
   3.7193718671916103`*^9, 3.7193718690751095`*^9}, 3.71937292455938*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"J", "=", 
   RowBox[{"GN", ".", 
    RowBox[{"Transpose", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "1"}], ",", "0", ",", 
         RowBox[{"-", "0.5"}], ",", 
         RowBox[{"-", "1"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", 
         RowBox[{"-", "2"}], ",", 
         RowBox[{"-", "2"}]}], "}"}]}], "}"}], "]"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.7193718198330717`*^9, 3.71937185386939*^9}, 
   3.7193720311259985`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"BB", "=", 
   RowBox[{
    RowBox[{"Inverse", "[", "J", "]"}], ".", "GN"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.7193725888799415`*^9, 3.7193726041348104`*^9}, {
   3.7193729403409953`*^9, 3.719372941875224*^9}, 3.7193739292149444`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"B", "=", 
   RowBox[{"{", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", "0", ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "2"}], "]"}], "]"}], ",", "0", ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "3"}], "]"}], "]"}], ",", "0", ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "4"}], "]"}], "]"}], ",", "0"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"2", ",", "1"}], "]"}], "]"}], ",", "0", ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"2", ",", "2"}], "]"}], "]"}], ",", "0", ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"2", ",", "3"}], "]"}], "]"}], ",", "0", ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"2", ",", "4"}], "]"}], "]"}]}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"2", ",", "1"}], "]"}], "]"}], ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"2", ",", "2"}], "]"}], "]"}], ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "2"}], "]"}], "]"}], ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"2", ",", "3"}], "]"}], "]"}], ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "3"}], "]"}], "]"}], ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"2", ",", "4"}], "]"}], "]"}], ",", 
       RowBox[{"BB", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "}"}]}], "}"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.7193727004907093`*^9, 3.7193727038266706`*^9}, {
  3.7193729545852346`*^9, 3.7193730506857367`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"DD", ":=", 
   RowBox[{
    FractionBox["EE", 
     RowBox[{"1", "-", 
      SuperscriptBox["v", "2"]}]], 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"1", ",", "v", ",", "0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"v", ",", "1", ",", "0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0", ",", 
        FractionBox[
         RowBox[{"1", "-", "v"}], "2"]}], "}"}]}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"EE", "=", 
   RowBox[{"3.", "*", 
    SuperscriptBox["10", "7"]}]}], ";", 
  RowBox[{"v", "=", "0.3"}], ";"}]}], "Input",
 CellChangeTimes->{{3.7193734249718533`*^9, 3.7193734815540285`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"ke", "[", 
    RowBox[{"eta_", ",", "psi_"}], "]"}], "=", 
   RowBox[{"Simplify", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Transpose", "[", "B", "]"}], ".", "DD", ".", "B"}], "/", 
      SuperscriptBox["10", "7"]}], "*", 
     RowBox[{"Abs", "@", 
      RowBox[{"Det", "[", "J", "]"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.719373411680382*^9, 3.719373417160321*^9}, {
  3.7193734838579955`*^9, 3.719373523692178*^9}, {3.7193737088921824`*^9, 
  3.7193737421085043`*^9}, {3.719373956416152*^9, 3.719373999857562*^9}, {
  3.7193744802070665`*^9, 3.7193744862783623`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", "@", 
  RowBox[{"UpperTriangularize", "@", 
   RowBox[{"ke", "[", 
    RowBox[{
     FractionBox[
      RowBox[{"-", "1"}], 
      SqrtBox["3"]], ",", 
     FractionBox[
      RowBox[{"-", "1"}], 
      SqrtBox["3"]]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.7193739813107395`*^9, 3.7193740053305473`*^9}, {
  3.719456523814316*^9, 3.7194565259203873`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1.2467344367125135`", 
      RowBox[{"-", "0.37258730292557873`"}], 
      RowBox[{"-", "1.1329832153767232`"}], "0.17758520920708112`", 
      RowBox[{"-", "0.33406148549318865`"}], "0.09983446692899887`", 
      "0.22031026415739838`", "0.09516762678949871`"},
     {"0.`", "0.6878534823241454`", "0.1450848602539982`", 
      RowBox[{"-", "0.36284999279331603`"}], "0.09983446692899893`", 
      RowBox[{"-", "0.18430978509969032`"}], "0.1276679757425817`", 
      RowBox[{"-", "0.14069370443113915`"}]},
     {"0.`", "0.`", "1.148222989324213`", "0.04991723346449944`", 
      "0.3035819375982098`", 
      RowBox[{"-", "0.0388753711390411`"}], 
      RowBox[{"-", "0.3188217115456993`"}], 
      RowBox[{"-", "0.15612672257945653`"}]},
     {"0.`", "0.`", "0.`", "0.4063922040718573`", 
      RowBox[{"-", "0.04758381339474937`"}], "0.09722536254260773`", 
      RowBox[{"-", "0.17991862927683125`"}], 
      RowBox[{"-", "0.14076757382114904`"}]},
     {"0.`", "0.`", "0.`", "0.`", "0.089511505260241`", 
      RowBox[{"-", "0.026750564790416875`"}], 
      RowBox[{"-", "0.05903195736526202`"}], 
      RowBox[{"-", "0.025500088743832666`"}]},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.04938565807461573`", 
      RowBox[{"-", "0.03420853099954093`"}], "0.037698764482466875`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.15754340475356302`", 
      "0.08645918453379046`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.2437625137698213`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.719373992821232*^9, 3.719374005948048*^9}, 
   3.7194565343930845`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"a", "=", 
   FractionBox["1", 
    SqrtBox["3"]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"kefinal", " ", "=", " ", 
   RowBox[{
    RowBox[{"ke", "[", 
     RowBox[{
      RowBox[{"-", "a"}], ",", 
      RowBox[{"-", "a"}]}], "]"}], "+", 
    RowBox[{"ke", "[", 
     RowBox[{
      RowBox[{"-", "a"}], ",", "a"}], "]"}], "+", 
    RowBox[{"ke", "[", 
     RowBox[{"a", ",", 
      RowBox[{"-", "a"}]}], "]"}], "+", 
    RowBox[{"ke", "[", 
     RowBox[{"a", ",", "a"}], "]"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.719374022040625*^9, 3.7193740619412537`*^9}, {
  3.719374366490471*^9, 3.7193743711331587`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", "[", 
  RowBox[{"UpperTriangularize", "[", "kefinal", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.719374810211156*^9, 3.719374816273259*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"2.7467244294167363`", 
      RowBox[{"-", "0.7417582417582416`"}], 
      RowBox[{"-", "2.458262890955198`"}], "0.24725274725274718`", 
      RowBox[{"-", "1.6768808114961957`"}], "0.6593406593406591`", 
      "1.3884192730346574`", 
      RowBox[{"-", "0.1648351648351648`"}]},
     {"0.`", "1.4898562975486047`", "0.16483516483516478`", 
      RowBox[{"-", "0.6656804733727809`"}], "0.6593406593406592`", 
      RowBox[{"-", "0.9763313609467456`"}], 
      RowBox[{"-", "0.08241758241758237`"}], "0.15215553677092136`"},
     {"0.`", "0.`", "2.6024936601859676`", "0.3296703296703295`", 
      "1.3884192730346574`", 
      RowBox[{"-", "0.08241758241758236`"}], 
      RowBox[{"-", "1.532650042265427`"}], 
      RowBox[{"-", "0.41208791208791207`"}]},
     {"0.`", "0.`", "0.`", "1.0777683854606928`", 
      RowBox[{"-", "0.1648351648351648`"}], "0.15215553677092142`", 
      RowBox[{"-", "0.41208791208791207`"}], 
      RowBox[{"-", "0.5642434488588336`"}]},
     {"0.`", "0.`", "0.`", "0.`", "3.8165680473372774`", 
      RowBox[{"-", "0.8241758241758239`"}], 
      RowBox[{"-", "3.528106508875739`"}], "0.32967032967032955`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "2.0033812341504644`", 
      "0.24725274725274718`", 
      RowBox[{"-", "1.1792054099746407`"}]},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "3.6723372781065082`", 
      "0.2472527472527472`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "1.5912933220625525`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.719374817030927*^9, 3.719456543397014*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"R", "=", 
   RowBox[{"{", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"i", "[", "0", "]"}], ",", 
       RowBox[{"i", "[", "1", "]"}], ",", 
       RowBox[{"i", "[", "2", "]"}], ",", "0", ",", "0", ",", "0", ",", "0", 
       ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"j", "[", "0", "]"}], ",", 
       RowBox[{"j", "[", "1", "]"}], ",", 
       RowBox[{"j", "[", "2", "]"}], ",", "0", ",", "0", ",", "0", ",", "0", 
       ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", 
       RowBox[{"i", "[", "0", "]"}], ",", 
       RowBox[{"i", "[", "1", "]"}], ",", 
       RowBox[{"i", "[", "2", "]"}], ",", "0", ",", "0", ",", "0", ",", "0", 
       ",", "0", ",", "0"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", 
       RowBox[{"j", "[", "0", "]"}], ",", 
       RowBox[{"j", "[", "1", "]"}], ",", 
       RowBox[{"j", "[", "2", "]"}], ",", "0", ",", "0", ",", "0", ",", "0", 
       ",", "0", ",", "0"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
       RowBox[{"i", "[", "0", "]"}], ",", 
       RowBox[{"i", "[", "1", "]"}], ",", 
       RowBox[{"i", "[", "2", "]"}], ",", "0", ",", "0", ",", "0"}], "}"}], 
     ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
       RowBox[{"j", "[", "0", "]"}], ",", 
       RowBox[{"j", "[", "1", "]"}], ",", 
       RowBox[{"j", "[", "2", "]"}], ",", "0", ",", "0", ",", "0"}], "}"}], 
     ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
      "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
       "0", ",", "0", ",", 
       RowBox[{"i", "[", "0", "]"}], ",", 
       RowBox[{"i", "[", "1", "]"}], ",", 
       RowBox[{"i", "[", "2", "]"}]}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
      "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
       "0", ",", "0", ",", 
       RowBox[{"j", "[", "0", "]"}], ",", 
       RowBox[{"j", "[", "1", "]"}], ",", 
       RowBox[{"j", "[", "2", "]"}]}], "}"}]}], "\[IndentingNewLine]", 
    "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"i", "[", "0", "]"}], "=", "0"}], ";", 
  RowBox[{
   RowBox[{"i", "[", "1", "]"}], "=", 
   RowBox[{"-", "1"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"i", "[", "2", "]"}], "=", "0"}], ";", 
  RowBox[{
   RowBox[{"j", "[", "0", "]"}], "=", 
   RowBox[{"-", "1"}]}], ";", 
  RowBox[{
   RowBox[{"j", "[", "1", "]"}], "=", "0"}], ";", 
  RowBox[{
   RowBox[{"j", "[", "2", "]"}], "=", "0"}], ";"}]}], "Input",
 CellChangeTimes->{
  3.7193743144168615`*^9, {3.7193743566155853`*^9, 3.719374356875285*^9}, {
   3.719374399667907*^9, 3.719374428308692*^9}, {3.719374633590417*^9, 
   3.719374636097826*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", "[", 
  RowBox[{"UpperTriangularize", "[", 
   RowBox[{
    RowBox[{"Transpose", "[", "R", "]"}], ".", "kefinal", ".", "R"}], "]"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.7193744381486397`*^9, 3.719374451016511*^9}, {
  3.7193745033107853`*^9, 3.7193745050608225`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1.4898562975486047`", 
      RowBox[{"-", "0.7417582417582416`"}], "0.`", 
      RowBox[{"-", "0.6656804733727809`"}], "0.16483516483516478`", "0.`", 
      RowBox[{"-", "0.9763313609467456`"}], "0.6593406593406592`", "0.`", 
      "0.15215553677092136`", 
      RowBox[{"-", "0.08241758241758237`"}], "0.`"},
     {"0.`", "2.7467244294167363`", "0.`", "0.24725274725274718`", 
      RowBox[{"-", "2.458262890955198`"}], "0.`", "0.6593406593406591`", 
      RowBox[{"-", "1.6768808114961957`"}], "0.`", 
      RowBox[{"-", "0.1648351648351648`"}], "1.3884192730346574`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "1.0777683854606928`", "0.3296703296703295`", 
      "0.`", "0.15215553677092142`", 
      RowBox[{"-", "0.1648351648351648`"}], "0.`", 
      RowBox[{"-", "0.5642434488588336`"}], 
      RowBox[{"-", "0.41208791208791207`"}], "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "2.6024936601859676`", "0.`", 
      RowBox[{"-", "0.08241758241758236`"}], "1.3884192730346574`", "0.`", 
      RowBox[{"-", "0.41208791208791207`"}], 
      RowBox[{"-", "1.532650042265427`"}], "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "2.0033812341504644`", 
      RowBox[{"-", "0.8241758241758239`"}], "0.`", 
      RowBox[{"-", "1.1792054099746407`"}], "0.24725274725274718`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "3.8165680473372774`", 
      "0.`", "0.32967032967032955`", 
      RowBox[{"-", "3.528106508875739`"}], "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "1.5912933220625525`", "0.2472527472527472`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "3.6723372781065082`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.7193744458308144`*^9, 3.719374451394889*^9}, {
   3.71937449766947*^9, 3.7193745054038687`*^9}, 3.7193746393026733`*^9, 
   3.719456552673996*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", "[", 
  RowBox[{
   RowBox[{"Simplify", "[", 
    RowBox[{"DD", ".", 
     RowBox[{"Evaluate", "[", 
      RowBox[{"B", "/.", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"eta", "\[Rule]", 
          FractionBox["1", 
           SqrtBox["3"]]}], ",", 
         RowBox[{"psi", "\[Rule]", 
          FractionBox["1", 
           SqrtBox["3"]]}]}], "}"}]}], "]"}]}], "]"}], ".", "R", ".", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
     RowBox[{
      RowBox[{"-", "1.17777"}], "*", 
      RowBox[{"10", "^", 
       RowBox[{"-", "6"}]}]}], ",", 
     RowBox[{
      RowBox[{"-", "9.66972"}], "*", 
      RowBox[{"10", "^", 
       RowBox[{"-", "6"}]}]}], ",", "0", ",", 
     RowBox[{"2.67425", "*", 
      RowBox[{"10", "^", 
       RowBox[{"-", "6"}]}]}], ",", 
     RowBox[{
      RowBox[{"-", "9.93532"}], "*", 
      RowBox[{"10", "^", 
       RowBox[{"-", "6"}]}]}], ",", "0"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7194565835937543`*^9, 3.7194566098513823`*^9}, {
  3.7194566583532524`*^9, 3.719456863196898*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {
       RowBox[{"-", "22.981850944270718`"}]},
      {
       RowBox[{"-", "42.02054844530335`"}]},
      {"2.5528959728400658`"}
     },
     GridBoxAlignment->{
      "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
        "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.5599999999999999]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}],
    Column], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.719456585015293*^9, 3.71945659582557*^9}, 
   3.719456723231927*^9, {3.7194567547409315`*^9, 3.719456772848966*^9}, {
   3.7194568217034636`*^9, 3.719456863867663*^9}}]
}, Open  ]]
},
WindowSize->{1776, 592},
WindowMargins->{{Automatic, 62}, {Automatic, 124}},
Magnification:>1.6 Inherited,
FrontEndVersion->"11.0 for Microsoft Windows (64-bit) (2016\:5e7410\:67088\
\:65e5)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 227, 4, 74, "Input"],
Cell[788, 26, 734, 21, 48, "Input"],
Cell[1525, 49, 559, 17, 74, "Input"],
Cell[2087, 68, 275, 6, 74, "Input"],
Cell[2365, 76, 2160, 63, 139, "Input"],
Cell[4528, 141, 710, 23, 110, "Input"],
Cell[5241, 166, 654, 16, 81, "Input"],
Cell[CellGroupData[{
Cell[5920, 186, 398, 12, 110, "Input"],
Cell[6321, 200, 2123, 44, 246, "Output"]
}, Open  ]],
Cell[8459, 247, 668, 21, 120, "Input"],
Cell[CellGroupData[{
Cell[9152, 272, 180, 3, 74, "Input"],
Cell[9335, 277, 2082, 44, 246, "Output"]
}, Open  ]],
Cell[11432, 324, 3167, 77, 469, "Input"],
Cell[CellGroupData[{
Cell[14624, 405, 310, 7, 74, "Input"],
Cell[14937, 414, 2810, 55, 350, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17784, 474, 1125, 33, 164, "Input"],
Cell[18912, 509, 930, 25, 148, "Output"]
}, Open  ]]
}
]
*)

